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Abstract. For particles emerging from a second order QCD phase transition, we show that a 
recently introduced shape parameter of the Bose-Einstein correlation function, the Levy index of 
stability equals to the correlation exponent - one of the critical exponents that characterize the 
behavior of the matter in the vicinity of the second order phase transition point. Hence the shape of 
the Bose-Einstein / HBT correlation functions, when measured as a function of bombarding energy 
and centrality in various heavy ion reactions, can be utilized to locate experimentally the second 
order phase transition and the critical end point of the first order phase transition line in QCD. 


INTRODUCTION 

The study fractal phenomena was initiated in high energy particle and nuclear physics 
by Bialas and Peschanski in ref. [1], with the motivation of searching for a second order 
phase transition by studying intermittency or the the power-law behavior of moments of 
the multiplicity distribution in narrowing bins of the momentum space, see refs. [2, 3] 
for excellent reviews on this topic. 

The mathematical properties of Bose-Einstein correlation functions for Levy stable 
sources were written up by three of us in refs. [4, 5], and are recapitulated in the next 
section. In ref. [6] we have added a physical interpretation and showed, that in case of 
jet physics, the fractal properties of QCD cascades can naturally be measured by the 
Levy index of stability of the Bose-Einstein correlation functions. Our analytic results 
were similar in spirit to the numerical investigations of Wilk and collaborators in ref. [7]. 
Note that these correlations are frequently referred to as Hanbury Brown - Twiss or HBT 
correlations in the literature of heavy ion physics. 

Bialas realized, that Bose-Einstein correlations and intermittency might be deeply 
connected [8], and considered a distribution of Gaussians where the radius parameter 
of the Gaussian has a power-law distribution, thus giving a way to the study fractals in 
coordinate space with the help of Bose-Einstein correlations. Brax and Peschanski were 
the first to introduce Levy distributions, in momentum space, to multiparticle production 
in high energy physics [9]. They have suggested to use the measured value of the Levy 
index of stability to signal quark gluon plasma production in heavy ion physics. Here we 
reconnect these seemingly different topics, and show how the excitation function of the 
shape parameter of the correlation function can be utilized to locate experimentally the 
critical end-point of QCD. 


BOSE-EINSTEIN CORRELATIONS & LEVY STABLE SOURCES 


The two-particle Bose-Einstein correlation function is defined with the help of the two- 
particle and single-particle invariant momentum distributions as: 


C2(ki,k2) 


V2(ki,k2) 
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( 1 ) 


If long-range correlations can be neglected or corrected for, and if the short-range cor¬ 
relations are dominated by Bose-Einstein correlations, this two-particle Bose-Einstein 
correlation function is related to the Eourier-transformed source distribution. Eor clarity, 
let us consider the case of a one-dimensional, factorized source. 


S{x,k)= f{x)g{k). (2) 

In this case [4, 5], the Bose-Einstein correlation function is 

C2{kuk2) = ^ + \hq)\\ (3) 

where the Eourier transformed source density (often referred to as the characteristic 
function) and the relative momentum are defined as 


f{q) 
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q = k\ —k 2 . 


(4) 


Eor the case of the jets decaying to jets to jets and so on, as well as at a second order 
phase transition, where fluctuations appear on all possible scales with a power-law tailed 
distribution, the final position of a particle is given by a large number of position shifts, 
hence the distribution of the final position x is obtained as a convolution, 

^ fix) = [ U^^idxilTj^Jjixj) 5(;c- Xk). (5) 

i=l k=l 

In the case of particle emission from QCD jets, that the fractal defining the particle 
emission is infrared stable: adding one more, very soft gluon does not change the result¬ 
ing source distributions. A similar property holds for systems at a second order phase 
transition: the system becomes invariant under a renormalization group transformation. 
Bose-Einstein correlation functions for such particle emitting sources were evaluated 
recently by three of us, which we summarize here following refs. [4] and [5]. 

Various forms of the Central Eimit Theorem state, that under certain conditions, the 
distribution of the sum of large number of random variables converges (for n oo) to a 
limit distribution. In case of “normal" elementary processes, that have finite means and 
variances, the limit distribution of their sum is a Gaussian. In case of random motion in 
a thermal medium, such position distribution corresponds to normal diffusion. However, 
near a second order phase transition point, fluctuations appear on all scales and the 
variance of the elementary process diverges, corresponding to the so-called anomalous 
diffusion. In this case, the system still may be invariant under convolution, and the shape 
of the limit distribution becomes independent from the number of elementary steps. 



f(s) Levy sources 




Figure 1. (left) Source functions for univariate symmetric Levy laws, as a function of the 
dimensionless variable s = r/R, for various values of the Levy index of stability, a. (right) 
Bose-Einstein correlation (or HBT) correlation functions for univariate symmetric Levy laws, 
for a fixed scale parameter of 7? = 0.8 fm and various values of the Levy index of stability, a. 


Stable distributions are preeisely those limit distributions that ean oeeur in General¬ 
ized Central Limit theorems. Their study was begun by the mathematieian Paul Levy 
in the 1920’s. The stable distributions can be given in terms of their characteristic func¬ 
tions, as the Fourier transform of a convolution is a product of the Fourier-transforms, 

/( 9 ') =n/'(^), ( 6 ) 

i=l 

and limit distributions appear when the convolution of one more elementary process 
does not change the shape of the limit distribution, but it results only in a modification 
of the parameters of the limit distribution. The characteristic function of univariate and 
symmetric stable distributions is 

/(<?)= exp (f^6-|y< 5 r|“), (7) 

where the support of the density function f{x) is (— 00 , 00 ). Deep mathematical results 
imply that the index of stability, a, satisfies the inequality 0 < a < 2, so that the 

source distribution be always positive. These Levy distributions are indeed stable under 

convolutions, in the sense of the following relations: 

fi{q) = exp{iqdi-\yiq\^), = exp (/^5- ly-?!”), (8) 

i=l 
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Thus the Bose-Einstein correlation functions for uni-variate, symmetric stable distri¬ 
butions (after a core-halo correction, and a re-scaling) read as 

C(^;a) = 1-|-?iexp(—. (10) 

Refs. [4] and [5] discuss further examples and details and generalize these results to 
three dimensional, hydrodynamically expanding, core-halo type sources, as well as to 
three-particle correlations. 

















THE ANOMALOUS DIMENSION OE QCD JETS AND BE/HBT 


In QCD, jets emit jets that emit additional jets and so on. The resulting fractal structure 
of QCD jets was related to intermittency and power-law dependence of multiplicity 
moments on the bin-size in momentum space with the help of a beautiful geometric 
interpretation of the color dipole picture in refs. [10, 11, 12], and an infrared stable 
measure on the parton states, related to the hadronic multiplicity distribution. These 
ideas were developed further in refs. [13, 14, 15], [16] as well as in refs. [17, 18]. 



Figure 2. The phase-space of QCD jets in the (y,K) plane, where K = \og{kl). (a) The phase 
space available for a gluon emitted by a high energy qq system is a triangular region in the (y, k) 
plane, (b) If one gluon is emitted at (yi,Ki), the phase space for a second (softer) gluon is given 
by the area of this folded surface, (c) The total gluonic phase space can be described by this 
multifaceted fractal surface [10, 11, 12]. 

A high energy qq system radiates gluons according to the dipole formula 

an = —^aya(|), (11) 

hence the phase-space for the emission of a gluon is given by the relation 

\y\<]^\n{s/k]_), (12) 

which corresponds to the triangular region in a (y, InA:^) diagram as shown in Fig. 2 
(a). If two gluons are emitted, then the distribution of the hardest gluon is described by 
eq. (11). The distribution of the second, softer, gluon corresponds to two dipoles, the 
first is stretched between the quark and the first gluon, and the second between the first 
gluon and the anti-quark. The phase-space available for the second gluon corresponds 
to the folded surface in Fig. 2 (b), with the constraint 2 < 1 ’the first gluon is 

assumed to be the hardest one. This procedure can be generalized so that the emission 
of a third, still softer gluon corresponds to radiation from three color dipoles, with n 
gluons emitted already the emission of the n -f 1-th gluon is given by a chain of n -f 1 
dipoles. Thus, with many gluons, the gluonic phase space can be represented by a multi¬ 
faceted surface as illustrated in Fig. 2 (c). Each gluon adds a fold to the surface, which 
increases the phase-space for softer gluons. (Note, that in this process the recoils are 
neglected, as is normal in leading log approximation). Due to its iterative nature, the 
process generates a Koch-type fractal curve at the base-line. The length of this base-line 
of the partonic structure on Figure 2 (c) is proportional to the particle multiplicity. This 




curve is longer, when studied with higher resolution: it is a fraetal eurve, embedded into 
the four-dimensional energy-momentum spaee, eharaeterized by the fraetal dimension 

or one plus the anomalous dimension of QCD [10, 11, 12]. With the help of the Lund 
string fragmentation picture, this fractal in momentum spaee is mapped into a fractal in 
coordinate spaee, and the eonstant of eonversion is the hadronie string tension, K 1 
GeV/fm. This mapping does not ehange the fraetal properties of the eurve. 

A walk, where the length of the steps is given by a Levy distribution, and the direetion 
of the steps is random, eorresponds to a fraetal eurve, in physieal terms it ean be 
interpreted as the path of a test partiele performing a generalized Brownian motion. This 
motion is referred to as anomalous diffusion and the probability that the test partiele 
diffuses to distanees r greater than a eertain value of |5| is given by P{r > |5|) «« kl ”• 
This relation is valid for anomalous diffusion in any dimensions. Thus the Levy index 
of stability a is the fraetal dimension of the trajeetory of the eorresponding anomalous 
diffusion [19]. When we apply this result to QCD, there are two key eonsiderations. 

First, if gluon radiation is negleeted, the qq system hadronizes as a l-i-l dimensional 
hadronie string, whieh has no fraetal strueture. If the gluon emission is switched on, the 
emission of gluon n from one of the n dipoles corresponds to a step of an anomalous 
diffusion in the plane transverse to the given dipole. Henee the anomalous dimension of 
QCD equals to the Levy index of stability of this anomalous diffusion. 



(14) 


Seeond, data on Bose-Einstein e orrelations are often determined in terms of the invariant 
momentum differenee — a /—{pi — PiY- Bose-Einstein eorrelation funetions that 
depend on this invariant momentum difference ean be obtained within the framework of 
the so-ealled x-model. This model assumes a broad proper-time distribution, H{x) and 
very strong eorrelations between eoordinate and momentum in all direetions, x^/x 
p^/mqy Henee {x\ — X 2 ){pi — pi) '^2fnv’ details. In this ease, 

the Bose-Einstein eorrelation funetion measures the Eourier-transformed proper-time 
distribution H in the following, unusual manner: 





(15) 


where niy'^ stands for the (transverse) mass of the pair for (two)- or more jet events. 
Erom this relation it follows, that — 2aLj.y. Thus we find the following relationship 
between the running QCD eoupling eonstant and the exponent of an invariant relative 
momentum dependent Bose-Einstein eorrelation funetion aBEc^ 

= ^a^c- ( 16 ) 

In ref. [6] we have eompared this leading log result to NA22 and UAl eorrelation data 
of refs. [22],[23] and found a reasonable agreement with these data. 



BOSE-EINSTEIN CORRELATIONS AT A SECOND ORDER QCD 

PHASE TRANSITION 


The main motivation behind the experimental and theoretieal program of high energy 
physies is to study the phase diagram of hot and dense hadronie matter. Aeeording to re- 
eent lattiee QCD ealeulations at finite temperature and baryon density, there exist a line 
of first order phase transitions that separates the hadronie and the quark-gluon plasma 
(QGP) state. This line of the first order phase transitions ends at the eritieal end-point 
(CEP), where the transition from hadron gas to QGP beeomes a seeond order phase tran¬ 
sition. Reeent lattiee QCD ealeulations loeated [24] this CEP at = 162 ±2 MeV and 
= 360 ± 40 MeV. Below these baryoehemieal values, the transition from a hadron gas 
to a QGP beeomes a eross-over, and at vanishing net baryon density the eritieal temper¬ 
ature beeomes = 164 ± 2 MeV (the errors are statistieal only). In this ealeulation, the 
quark masses were already at the physieal value, but the eontinuum extrapolation was 
missing. S. Katz presented improvements at the Quark Matter 2005 eonferenee [25], 
using physieal quark masses and working towards the eontinuum extrapolation. He re¬ 
ported Tc— 189 ± 8 MeV for the eritieal temperature at hb = 0. 

At the CEP, the seeond order phase transition is eharaeterized by the fixed point 
of the renormalization group transformations. In a quark-gluon plasma, the vaeuum 
expeetation value of the quark eondensate c = (^q) vanishes , while in the hadronie 
phase, this vaeuum expeetation value beeomes non-zero. The eorrelation funetion of 
the order parameter is defined as p(R) = (c(r-t-R)c(r)) — (c)^ and measures the spatial 
eorrelation between the pions. At the CEP, this eorrelation funetion deeays as 

p(R) (17) 

a power-law. The parameter p is ealled as the exponent of the eorrelation funetion. 

Eor Eevy stable sourees, eorresponding to an anomalous diffusion with large fluetua- 
tions in eoordinate spaee, the eorrelation between the initial and aetual positions deeays 
also as a power-law, where the exponent is given by the Eevy index of stability a as 

p(R) (18) 

As we are eonsidering a QCD phase transition in a J = 3 three-dimensional eoordinate 
spaee we find that the eorrelation exponent equals to the Eevy index of stability, a = Tj. 
Stephanov, Rajagopal and Shuryak pointed out [26], that the universality elass of the 
seeond order QCD phase transition is that of the 3d Ising model. Eor this universality 
elass, the eorrelation exponent has been determined by Rieger [27] as 

a(Eevy) =ri (3d Ising) =0.50 ±0.05. (19) 

Eig. I indieates that the ehange in the shape of the eorrelation funetion is rather signif- 
ieant, if a deereases from its Gaussian value of 2 to 0.5, its eharaeteristie value at the 
2nd order QCD phase transition point. Eig. 3 illustrates how this shape parameter of 
the Bose-Einstein eorrelation funetion may depend on the relative temperature near the 
eritieal point. Henee studying the shape parameter of the two-partiele Bose-Einstein or 
HBT eorrelation funetions as a funetion of the bombarding energy or the eentrality of 



the heavy ion eollisions, a previously unknown tool is obtained to determine if the pions 
are emitted from the neighborhood of the eritieal end point of the QCD phase diagram. 

Furthermore, based on an universality elass argument, we have determined that the 
seeond order QCD phase transition at the eritieal end point will be signaled with the 
value of a = 0.5, a very spiky Bose-Einstein eorrelation funetion indeed. 


Levy index of stability 



Figure 3. Illustration of the behavior of the Levy index of stability of Bose-Einstein correlations 
as a function of the dimensionless temperature variable X = {T — Tc)/Tc in the neighborhood of 
the critical endpoint of the 1st order phase transition line in QCD. At the critical endpoint, the 
phase transition becomes 2nd order and the Levy index of stability decreases to the correlation 
exponent of QCD. As this transition has the same universality class as that of the 3d Ising 
model, one expects a decrease from the a 2 values that are characteristic to a Boltzmann 
gas and normal diffusion to a = 0.5, corresponding to the correlation exponent of QCD at the 
critical endpoint. As shown in Fig. 1, such a change in the shape parameter makes the Bose- 
Einstein correlation functions much sharper than a simple Gaussian, so the spiky structure of the 
correlation function could be used to search for this point experimentally. 


CONCLUSIONS 

We have reeapitulated earlier results that indieate, that the general shape of the Bose- 
Einstein or HBT eorrelation funetions is a stretehed exponential or Eevy stable form, 
where the Eevy index of stability beeomes a new shape parameter of the eorrelation 
funetion with 0 < a < 2 and the popular Gaussian parameterization eorresponds to the 
a = 2 partieular, speeial ease. Then we have studied two physieally interesting examples. 

In ease of partiele emission from jets, we have reeapitulated the eonneetion between 
the stability index of the Bose-Einstein/HBT eorrelation funetions and the running 
eoupling eonstant of QCD. 

We have also eonsidered a seenario, when the power-law tail of a Eevy distribution of 
the partiele emission in the coordinate space appears due to a seeond-order QCD phase 





transition. In this case, the Levy index of stability of the Bose-Einstein or HBT eorre- 
lation funetion was shown to be equal to the eorrelation exponent of QCD. This value 
is known to be 0.5 ± 0.05 from universality elass eonsiderations. Henee by measuring 
the exeitation funetion of the Levy index of stability (the shape parameter of the two- 
partiele Bose-Einstein or HBT eorrelation funetions), one ean experimentally determine 
the bombarding energy and eentrality range where a heavy ion eollision hits the eritical 
end point of QCD. Clearly, more work is neeessary to eheek to what extent this interest¬ 
ing effeet ean be masked by the deeays of various resonanees, hydrodynamie expansion, 
and by the time evolution of the partiele emitting souree between the seeond order phase 
transition point and the freeze-out temperature. 
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